Superluminal Neutrinos at OPERA Confront Pion Decay Kinematics 
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Violation of Lorentz invariance (VLI) has been suggested as an explanation of the superluminal 
velocities of muon neutrinos reported by OPERA. In this note we show that the amount of VLI 
required to explain this result poses severe difficulties with the kinematics of the pion decay, ex- 
tending its lifetime and reducing the momentum carried away by the neutrinos. We show that the 
OPERA experiment limits a = (v v — c)/c < 4 x 10 -6 . We then take recourse to cosmic ray data 
on the spectrum of muons and neutrinos generated in the earth's atmosphere to provide a stronger 
bound on VLI: (v -c)/c< 10~ 12 . 
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The recent OPERA report of superluminal veloc- 
ities for the muon neutrinos, v(v^)/c = 1 + a = 
2.5 x 10 -5 , has evoked much interest. Indeed present in- 
formation on neutrino oscillations suggests much stronger 
bounds on putative superluminal anomalies for neutrinos 
[3, El] • Still this recent experiment andprevious measure- 
ments at Fermilab Q and MINOS [EJ supporting this 
result prompted many theoretical and phenomenologi- 
cal comments. These possibilities include speculations of 
segregating the effect only into the za, sector @-@|. In 
this paper we study the implications of the superluminal 
velocities of the neutrinos on the kinematics of pion decay 
and show that superluminal velocities for are severely 
constrained by these considerations. The constraints de- 
rived here are not restricted to any specific model but 
merely probe into consequences of superluminal motion 
of from pion, kaon and other decays. 

Most of the attempts to understand the OPERA re- 
sult consider violation of Lorentz invariance (VLI) at the 
phcnomcnological level There are also theo- 

retical motivations stemming from string theory and from 
models with extra dimensions. In these models VLI in- 
creases with energy as a power law and have the general 
characteristic of modifying the maximum attainable ve- 
locity of the particles. 

The phenomen ology o f these models have been exten- 
sively studied 0, lldt- 12[. and important constraints on 
the level of VLI have been established. Of particular in- 
terest is the work of Cohen and Glashow [l3| , who discuss 
the possibility of — \ e~ e + or — \ -\-v G -\- u e 
and derive strong constraints on VLI. Other ideas of 
emission of gravitational radiation have also been dis- 



cussed [I4j. Keeping these in mind, additional assump- 



tions are required to accommodate the large superlumi- 
nal velocities reported by the OPERA collaboration. The 
very severe constraints come from the neutrino sector: 
neutrino-oscillation experiments suggest that the amount 
of Lorentz non-invariant contribution for all the three 



neutrinos to be the same (a Ve = = ot VT = a v ), 
as noted by Coleman and Glashow [J H|. The obser- 
vations of neutrinos from the supernova SN1987A [L5l | 
require that \a v \ < 10 -9 . The recent OPERA claim of 
a Vfi = 2.5 x 10 -5 , together with the supernova SN1987A 
constraints seems to indicate that the VLI parameter 
grows rapidly with energy, as suggested by some mod- 
els. 

In this paper wc note that such a large value of 
cv^, whether energy-independent or energy-dependent, 
will have many other phcnomcnological manifestations. 
Specifically, they would affect the kinematics and the 



rate of n — > /i 



decay, for high energy pions in ways 



that many experiments (OPERA included) would have 
detected. Moreover, the change in the rate of pion de- 
cay would affect the flux of the cosmic ray muons and 
muon neutrinos significantly, in conflict with observa- 
tions which extend up to ~ 10 4 GcV and ~ 10 5 GeV, 
respectively. In the present analysis we assume that the 
neutrinos detected at Gran Sasso arise exclusively from 
pion decay, even though there could be some contribu- 
tions from kaons. As a justification of this assumption we 
note that the charged kaon multiplicity at these energies 
is low ^0.3, much smaller than the pion multiplicity of 
~ 6, and that the transverse momenta of kaons are larger 
than that of the pions and the transverse momenta of 
neutrinos arising in X-dccay are larger than those from 
7r-dccay, so that the probability of K — > fi + con- 
tributing to the detector at Gran Sasso 730 kms away 
is reduced. Furthermore, the considerations presented 
here are equally applicable (with some numerical modi- 
fications) to -ftT-decay as well. A more detailed analysis 
of the kaon contributions is certainly warranted both in 
the context of OPERA results and for the analysis of 
cosmic-ray data, but will not change the conclusions of 
this paper significantly. 

In the formalism for VLI, given by Coleman and 
Glashow 0] , different particles achieve different terminal 



2 



velocities, and accordingly, for the discussion of 7r-decay, 
we make the minimal assumption that muon neutrinos 
have superluminal motion and the n , /x ± , being charged 
particles, have terminal velocities equal to the velocity 
of light to avoid Cerenkov radiation in vacuum. Thus 
unlike the analysis reviewed in the introduction, our 
analysis presented here does not directly apply to v e and 
v T , except indirectly because of neutrino oscillations. 
However, the two-body kinematics presented here, with 
appropriately chosen a, is valid in all cases where one 
of the emergent particles has a superluminal terminal 
velocity. In models, where a increases with energy, the 
constraints derived in this paper become much more 
stringent. Wc make the following standard assumptions: 

Al Energy-momentum conservation. 

A2 The relation dE/dp = v between the velocity 
of a particle and the change of its energy with 
momentum. This classical relation applies also to 
the group velocity of waves, v group = duj/dk, and 
extends it to wave mechanics as well. 

A3 The positivity of energy for free particles (by which 
we exclude Tachyons). The assumed E—p relations 
for different fields are variants of deformed forward 
light cones or mass hyperboloids. 

These criteria are applicable to most existing VLI models. 
The assumption A2 for the muon neutrinos implies 



dE 



v(p)dp, 



(1) 



where v(p) > 1 beyond some small value of momentum 
Pmin that is much larger than the tiny sub eV mass of the 
neutrinos, m v , which we neglect. Defining, in general, 



/8E\ , 



(2) 



as the effective average over the neutrino momenta de- 
tected in the OPERA experiment, centered around 17 
GeV, we write the energy-momentum relation at high 
energies as 



E v =p u (l + a) , 



(3) 



where a corresponds to VLI, required to understand the 
OPERA anomaly. The kinematic analysis begins with 
the standard mass-energy relation for n, /j: 

25, = (p? + m*f* . (4) 

We then express the four vector of the decaying pion as 

(K,P*,0,0) (5) 

and those of the final neutrino and muon as 

{E v ,p vl ,p vt ,0) and (E^p^p^t, 0) , (6) 



respectively, where the longitudinal components of mo- 
menta are taken to be p v i = rjp^ and p^e = (1 — r/)p w and 
the transverse components as p vt = —p^t = Pt ■ With this 
choice the conservation of all the spatial components of 
momenta is evident. 

We still need to satisfy the energy conservation: 



E-n — E v + E,j, , 



(7) 



with: 



K = [pi + ml] 1/2 ; E v 



[pW + Pt} 1 ' 2 (1 + ") 



and 



E„ = [pl(l-ri) 2 +P 2 t+ml] 1/2 



(8) 



Keeping in mind that in accelerator experiments in- 
cluding OPERA, m-x/pir, rn^/pn and pt/p-K are very 
small, wc expand the square root and keep only the lead- 
ing term to get 



2p w 



P\ 



m. 



2p,r(l-T?) 

Rearranging we can write 



P 2 (l + «) 
Zp-kV 



Pt 



pi 



2pl 2p2(l_ ?y ) 2plr, 



Spiv 



(9) 



(10) 



1 



or a 



2pW+pi 



2 2 V 



Pt 



(1-7/) (1 - 77) 



Since p\ is positive, this yields a constraint: 

™2 



1 

2p~l 



(11) 



(12) 



In the OPERA experiment the typical energy of the neu- 
trinos is ~ 17 GeV that arise from the decay of pions 
with a mean energy of ~ 60 GeV, so that the typical 
value of (77) p» 0.3. Inserting this value of 77 into equation 
[T^we obtain the bound: 



a OPERA 



< 



1 



0.6p2 



0.7 



4 x 10" 



(13) 



Note that this bound on the superluminal parameter, a, 
is significantly smaller than 2.5 x 10 -5 estimated from 
the time profiles of the events and the GPS timing in 
their experiment. 

Next, we address the question, whether 77 could indeed 
be smaller than the assumed value of ~ 0.3 which would 
allow the value of a, estimated in our analysis, to be 
consistent with the OPERA result. For this, special se- 
lection effects should conspire to push 77 down to ~ 0.05. 
We note that this hypothesis would imply significant en- 
hancement of the lifetime of the pions. To see this, note 
that within this standard kinematics of pion decay, the 



3 



value of the rj parameter is uniformly distributed in the 
range 



< rj < 1 - 



(14) 



i.e., in the range ~ — 0.5. The phase space for the pion 
decay is directly proportional to this range and any re- 
duction in this range will have a corresponding reduction 
in the rate of decay. It is straightforward to perform the 
Lorentz non-invariant phase space integral after modi- 
fying the (5-functions representing mass shell conditions 
according to equations [3HS1 Such a calculation yields an 
integral directly proportional to rj max . Thus with the re- 
duction of (77) to 0.05, the pion lifetime will be extended 
by a factor of 6 or more, which is excluded by various 
accelerator experiments, including OPERA. 

As seen clearly from equations [TmiT^l the bounds get 
stronger in proportion to y\ , or even with higher powers 
in models with a increasing with energy invoked recently 
for explaining how the the OPERA results need not be 
flavor specific and still be consistent with the small a in- 
ferred from SN1987A neutrinos. Accordingly the bounds 
on VLI become extremely stringent for the ultra-high en- 
ergy muons and neutrinos observed in deep underground 
experiments at Kolar Gold Fields, Kamiokande, Baksan, 
Ice-Cube and other experiments [16l4l8| . 

Before we discuss these cosmic ray observations, we 
note that the fraction of the momentum carried away 
by the muon in the standard decay kinematics of the 
pion, (1 — rj) is in the range of <~ 0.5 — 1. The spectrum 
of muons generated in the Earth's atmosphere is well 
measured up to energies of ~ 10 5 GeV and we confine our 
analysis to the spectrum up to ~ 4 x 10 4 GeV where the 
muons arise mainly from the decay of pions and kaons 
and the contributions of muons generated by neutrino 
interactions in rock to the depth intensity curve could be 
neglected. The observed differential energy spectrum is 
well represented by the theoretical estimate [19f : 



A K (l 



1 - r))iz£n 
1 - T]) K £ K 



E+(1-7]) k £k 



E-e,(i5) 



where: 

/3 = spectral index of the cosmic ray spectrum ~ 2.65 
£ v = h o (0)m v i 'cTk 
£k = h (8)m K /cr K 

ho(9) = 7 x 10 5 sec# cm, the scale height of the Earth's 

atmosphere at a zenith angle 9 

t tt/k — decay lifetimes of pions/kaons at rest 

(1 — r))x/K = the fractional momenta carried away by the 

muons in pion/kaon decay averaged over the spectrum of 

cosmic rays, around the energy band of interest 

A v /k = Constants. 



These constants are estimated from the inclusive cross 
sections for the production of pions and kaons at high en- 
ergies and indicate that the net contribution of K-dec&y 
is ~ 10% for the muons and about ~ 70% of the total flux 
of neutrinos at the highest energies. A similar expression 
for the flux of neutrinos generated in the atmosphere re- 
sults when we replace (1 — 77) by 77 in eciuation ll5l Notice 
that at very high energies > 10 3 GeV, with E » £ v /k> 
the spectra of muons and neutrinos become steeper with 
a spectral index ~ ((3 + 1). Furthermore, the spectral 
intensities became proportional to (1 — ?])^ or {rf)P as the 
case might be. 

Now the spectrum of muons presented by Novoseltsev 
et al [13] , fits well with equation [TS1 that assumes that 
r n / K are constants. Thus (77) has to be constant up to 
energies of ~ 4 x 10 4 GeV. Note that equation [15] is sen- 
sitive to change in (77) in two ways - first through the 
change (1 — 77}^ and more importantly through its effect 
on extending the lifetime of pions and kaons. Thus the 
muon data imply 



ar) < 10 



-11 



(16) 



Much more extensive data of the atmospheric muons 
(2 x 10 10 events) and upward neutrinos (2 x 10 4 events) 
of energies in the range of 1 — 400 TeV, generated by 
energetic cosmic rays from the other side of the Earth 
have been provided recently by the south pole Ice-Cube 
experiment [Hj], which shows a good fit with an in- 
dex ~ (/3 + 1) ~ 3.65 at energies E » £ 7T /k- Thus their 
observations imply a constraint 



ar) < 10 



-13 



(17) 



Keeping in mind that we can not allow significant changes 
in rj as they will affect the spectral slope and spectral 
intensities of the muons and neutrinos, the limits derived 
here represent bounds on the supcrluminal parameter a, 
which may be stated conservatively as 



a < 10 



-12 



(18) 



allowing nearly a factor of ten variance for any contribu- 
tions of the uncertainties in the cosmic-ray observations 
and the approximations used in our analysis. It should 
be noted that since spectra of both the muons and the 
neutrinos are fit very well the estimates, which assume 
and 77 to be constants, the bound on the VLI param- 
eter a follows exclusively from kinematic considerations. 
Indeed accurate spectra of atmospheric muons and/or 
neutrinos even at lower energies of a several TeV can be 
used to improve the limits presented here. Our limit on 
a v is far more stringent than that derived from the ob- 
servations of the neutrinos emitted in SN1987A. It may 
be appropriate to note here that the observation of even 
a single event initiated by the GZK neutrinos in detec- 
tors like ANITA [2(J would improve the bound on a to 



4 



~ 1CU 20 . To assess the impact of these results on spe- 
cific models, we note that in general the matrix element 
in VLI theories may also have novel energy dependence, 
however they are unlikely to exactly cancel the above 
purely kinematic effects derived in this paper. 

We would like to draw attention to an independent 
analysis of the Ice-Cube data by Bi et al [2l[ , who assume 
that the superluminal may be treated as having an ef- 
fective mass, m e ff = [m 2 + 2ap 2 ] 1/ ' 2 , so that the decay 



mode 



/' 



becomes forbidden beyond a thresh- 



old energy for the neutrinos. This analysis yields result 
similar to our results, which we have derived showing the 
progressive kinematic restriction of the phase space avail- 
able for 7r-decay, leading to a monotonic increase of pion 
life-time with energy. 

In summary, we presented here a strong constraint 
(v — c)/c < 10~ 12 on the amount of violation of Lorentz 
invariance from pion decay kinematics and cosmic ray 
data. Careful observations of the fluxes of very high en- 
ergy muons and neutrinos at accelerators and in cosmic 
rays, and their comparison with the expected fluxes will 
constrain any possible variation of the decay life-time of 
the pion, which in turn, will lead to better bounds than 
those reported here. 
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